Pore formation in fluctuating membranes 



O 

o 

(N 

I— J 



Oded Farago^'U and Christian D. Santangelo^'O 

^Materials Research Laboratory, University of California, Santa Barbara, 9310^ 

^Department of Physics, University of California, Santa Barbara, 93106 

(Dated: February 2, 2008) 

We study the nucleation of a single pore in a fluctuating lipid membrane, specifically taking into 
account the membrane fluctuations, as well as the shape fluctuations of the pore. For large enough 
pores, the nucleation free energy is well-described by shifts in the effective membrane surface tension 
and the pore line tension. Using our framework, we derive the stability criteria for the various pore 
formation regimes. In addition to the well-known large-tension regime from the classical nucleation 
theory of pores, we also find a low-tension regime in which the effective line and surface tensions can 
change sign from their bare values. The latter scenario takes place at sufficiently high temperatures, 
where the opening of a stable pore of finite size is entropically favorable. 
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I. INTRODUCTION 

Lipid bilayers play an important role in living cells as 
barriers separating the inside of the cell from the extra- 
cellular environment, as well as segregating the cell into 
separate internal compartments [l|. A vital feature of 
those membranes is the ability to remain intact under 
variety of external perturbations 1^. At the same time, 
however, many cellular processes, including endo- and 
exocytosis, lysis, and cell signaling require breaking the 
membrane structure and forming a non-bilayer transient, 
such as a membrane pore 00, IS, 6] . The opening of sta- 
ble pores in biological membranes is also an important 
step for drug delivery Q and gene therapy [3|. Conse- 
quently, much attention has been focused on understand- 
ing the processes leading to the formation of pores and 
the mechanisms controlling their stability. 

Two types of pores dominate membrane permeability: 
free lipid pores and peptide-lined pores [9|. The interest 
in lipid pores has greatly increased in the past few years 
with the development of new experimental techniques to 
induce and study pore formation in biomimetic, single- 
component, lipid membranes. One approach to nucle- 
ate a pore is known as electroporation, where an elec- 
tric field that produces compressive stress is used to dis- 
rupt the membrane \^ Jd, J^, J_3, J^] . Holes in lipid 
membranes have also been opened by other methods, in- 
cluding intense illumination |15L I1b | , suction through a 
micropipet te l lTl . adhesion on porous or decorated sub- 
strates |iMll3' ^'^d osmotic swelling |20ll2ll |. 

Most theories of pore formation to date derive from 
a model based on classical nucleation theory [23 ■ The 
model conceives the membrane as a two dimensional elas- 
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tic medium characterized by a free energy per unit area 
( "surface tension" ) a. The formation of a circular hole of 
radius tq is driven by the reduction in the tension energy 
— cTTrrg and is opposed by an edge energy proportional 
the pore perimeter r27rro, where F, the line tension, de- 
notes the energy per unit length along the pore's rim. 
The net energy is, thus, given by 



E = T 27rro — cmrQ . 



(1) 



Assuming F > and cr > 0, Eq.QJ predicts that a pore 
with a radius larger than the critical value of 



F 
ro > - 



(2) 



is unstable in the sense that it will grow without bound 
and, ultimately, will rupture the membrane. Such a large 
pore will be created only if the nucleation energy barrier 



A£; = 



^F^ 



(3) 



is accessible by thermal fluctuations. For typical esti- 
mates of the fine tension, F ~ 10"*^ dyn [M I2I ^. 
thermally driven rupture requires a surface tension on 
the order of 1 dyn/cm. 

The problem with the above model is that it pre- 
cludes the existence of stable pores of finite size. Nev- 
ertheless, long-lived pores that remained open for sev- 
eral seconds before resealin g ha ve been observed in ex- 
periments HI El S m I23. Opening of transient 
pores has been also reported is several computer sim- 
ulations m 113, m lla]. This has led people to re- 
examine the basic assumptions underlying Eq.J^I- Im- 
proved theoretical models succeeded in explaining the 
formation of stable or long-lived metastable pores by 
considering the fact that, once pores have been nucle- 
ated, their further opening is expected to relax the sur- 
face tension j24, 28, 30, 31, 32, 33, 34] . In the case of a 
planar membrane, it is the increase of the area density 
of the lipids (occurring concurrently to the dilation of 



the pore) that reduces the mechanical tension. For vesi- 
cles, the opening of a pore allows the internal contents 
to escape, which reduces the osmotic pressure and the 
associated Laplace tension. 

Most of the theoretical models have so far, however, ne- 
glected the effect of membrane fluctuations on the open- 
ing and thermodynamic stability of pores. The entropy 
associated with the shape of the pore has mostly been 
ignored also. A few recent studies of the entropic contri- 
bution to the free energy of nucleating a pore, have led to 
some new interesting predictions. The most remarkable 
result has been obtained by Shillcock and Boal [35] in 
computer simulations of two-dimensional fluid tethered 
surfaces. They found that pores appeared at zero, and 
even small negative surface tension. Their interpretation 
of this surprising finding was that entropy, which favors 
the formation of non-circular pores, reduces the effective 
line tension of the pore and makes it negative at suffi- 
ciently high temperature. An entirely different fluctua- 
tion effect has been discussed by Sens and Safran jSg 
who considered circular pores, but allowed membrane 
fluctuations. Their study suggests that positive stress 
must be applied in order to facilitate the opening of a 
pore in a fluctuating membrane, and that the nucleation 
barrier for pore formation is too high to be overcome by 
thermal fluctuations. More recently, we [33 have demon- 
strated that the primary effect of membrane fluctuations 
on circular holes is to reduce the effective surface ten- 
sion, thereby making the opening of a pore more difficult 
in comparison to the zero-temperature case. 

In this paper, we carry out a statistical mechanical 
analysis of pore formation in bilayer membranes. Mem- 
brane elasticity is described by the Helfrich Hamiltonian, 
which includes the curvature energy and a surface ten- 
sion term. A line tension term is introduced to account 
for the energetic penalty at the pore edge. We calculate 
the free energy for nucleating a single pore, systemat- 
ically taking into account both membrane fluctuations 
and the entropy due to pore shape. We show that, for 
large enough pores, the pore nucleation free energy takes 
the form of Eq. (^ , with the surface and line tensions re- 
placed by their effective (renormalized) values, CTeff and 
Feff. The effective coefficients are usually smaller than 
the bare counterparts and, at high enough temperatures, 
may even be negative. When Fofr < 0, the opening of 
a stable pore becomes entropically favored and may oc- 
cur in weakly stressed membranes. The size of a thermal 
pore can be varied by changing the tension applied on the 
membrane, and the membrane is ruptured when ctcS > 0. 

The paper is organized as follows: The system Hamil- 
tonian is constructed in section ^1 We show that Eq.lQJ 
for the energy should augmented by terms represent- 
ing corrections due to the pore's shape and membrane's 
height fluctuations. Tracing over the relevant variables 
yields the corresponding free energy. The derivation of 
the free energy is presented in sections IIIII and IIVI in 
which the thermal corrections to the surface and line 
tensions are calculated. In section we discuss our re- 



sults and suggest some possible further extensions of the 
present study. 



II. DERIVATION OF THE HAMILTONIAN 

We consider a bilayer membrane consisting of N lipids 
that spans a planar circular frame of total area Ap = 
TrLp, in which a quasi-circular pore has been formed. For 
a nearly flat membrane with arbitrary parametrization 
X(2;l,a;2), the Helfrich Hamiltonian is given by 



n 



M 



dxidx2 y/g 



-H' 



kK 



(4) 



where g is the determinant of the metric tensor gaf} = 
da'K.-d/^'X-, while H and K denote the total and Gaussian 
curvatures, respectively. The elastic coefficients appear- 
ing in the Helfrich Hamiltonian are the surface tension a, 
the bending rigidity k, and the Gaussian rigidity R. We 
assume that the bilayer membrane is symmetric with no 
spontaneous curvature. The integration in Eq.Q is car- 
ried over the two-dimensional manifold M, representing 
the surface of the membrane. 



A. Gaussian curvature 

Understanding the contribution of the Gaussian cur- 
vature term [last term in Eq.l^J] to the free energy re- 
quires looking at the structure of the membrane on the 
molecular level. Two distinct models of pores have been 
discussed in the literature and are shown schematically in 
fig-D Fig.n(a) depicts a cylindrical pore where the lipids 
in the vicinity of the pore remain oriented parallel to the 
membrane surface. Such a pore is called "hydrophobic" , 
and the origin of the pore line tension is the energy due 
to the exposure of the tails of the lipids at pore's rim to 
water. The other case is of a "hydrophilic" pore, shown 
in fig.n(b), where the lipids curve at the pore's rim thus 
shielding their hydrophobic parts from the aqueous con- 
tact. The line tension of a hydrophilic pore is due to 
the curvature energy involved in the reorientation of end 
molecules. 

The Gaussian curvature term in the Helfrich Hamil- 
tonian is calculated differently for hydrophobic and hy- 
drophilic pores. For hydrophobic pores, this term re- 
duces to an integral of the geodesic curvature on the pore 
boundary jSg- A similar situation is encountered in the 
case of proteins and other membrane inclusions, where 
the orientation of the lipids at the membrane-inclusion 
boundary is determined by the structure of the protein 
[37[ . The change in the Gaussian curvature in that case 



ATYgauss = Sttk (cos V, - 1) , (5) 

where fl is the contact angle along the boundary. 



The case of a hydrophilic pore, on the other hand, is 
more subtle. Each of the two monolayers making up a 
bilayer can have an associated Helfrich energy, and the 
two monolayers are joined at the pore boundary. This, 
however, does not present a topological boundary and, 
therefore, there is no integral of the geodesic curvature. 
The opening of a hydrophilic pore results in a contribu- 
tion to the Hamiltonian due to the change in topology: 
The bilayer is topologically a sphere (if we assume both 
bilayers are linked at the outer radius Lp) , which changes 
genus upon the opening of a pore, becoming a torus. For 
a manifold without boundary, the Gauss-Bonnet theorem 
ensures that the total Gaussian curvature is a topological 
invariant and thus measures, to some extent, the global 
properties of the membrane. The change in Gaussian 
curvature energy due to the formation of a hydrophilic 
pore is given by 



An^ 



-Attk. 



(6) 





(a) 





(b) 



FIG. 1: Schematic microscopic models for hydrophobic (a) 
and hydrophilic (b) pores. For a hydrophilic pore, the bound- 
ary actually joins the two monolayer of the lipid membrane. 



The Gaussian curvature modulus may take both posi- 
tive or negative values, hence leading to either an increase 
(k < 0) or decrease {R > 0) in the energy upon the open- 
ing of a pore. Strict comparison with classical nucleation 
theory (Eq^ is possible only for k = and, therefore, 
we will restrict the following discussion to this special 
case. The Gaussian curvature term can be interpreted as 
an additional contribution that lowers or raises the free 
energy depending on the sign of k. This contribution is 
independent of the pore size. It will influence the proba- 
bility of opening a small nucleation pore, but will have no 
effect on the size to which (meta)stable long-lived pores 
grow. 



B. Membrane fluctuation energy 

The first two terms of the Helfrich Hamiltonian Q 
are local in character. Therefore, the above argument re- 
garding the absence of boundaries in a porous membrane 
fails, and the pore can be treated as if representing the 
inner membrane boundary. In order to study the statis- 
tical mechanical behavior of the membrane, we define a 
coordinate system {r,6) — (ro < r < Lp,0 < 9 < 2-k) in 
which the pore is described by a curve of constant r = ro : 



X(r,( 



[rcos(6') +?7a;(r, ( 
+h{r,e)z. 



]x + [rsin(0) 



(7) 



The function h{r,9) represents the height of the mem- 
brane above some flat reference plane. The function 
r cos{9)x + r sin(0)y -|- ff(r^ 9) is a mapping from coordi- 
nates (r, 9) in which the membrane pore will be circular, 
to points in three-dimensional space in which the pores 
will have an arbitrary shape (see fig. [SI- Thus, ff{ro,9) 
is a measure of the deviation of the pore from having a 
circular projected area, which we will assume to be small. 
Our choice or ro, which we will define as the radius of 
the quasi-circular pore, is made by equating the projected 
area of the pore to nrg . More specifically, we will require 
that ?7 satisfies the following boundary conditions (BCs): 



9-rf{ro,9) = 



27r 



d9 [ro + r • 77(ro,( 

On the outer (frame) boundary we set 
Tf{Lp,9)^0. 



(8) 
(9) 

(10) 



With the embedding defined by Eq.(|7I) we have, keep- 
ing terms up to quadratic order in ff and h, 

grr = I + 2f ■ drf] + [dr-qf + [drhf 

gre = f-d0ff+r9-drrj+drhd0h + drff-d0ff (11) 
gee = r^ + 2r9 ■ deif + {dgrf)^ + [dehf 

for the metric, thus giving us after a lengthy but straight- 
forward calculation 



y/g^r 



m 



W -rf+W ■T]2 



where 



m = iVxderiy) (f/r) - [rj^drriy) t 



(12) 



(13) 



One also finds that to lowest order the total curvature is 
given by 

H = -^da (Vgg'^^dpX) « z\/^h + 0{h\Tf). (14) 
v-9 




FIG. 2: The mapping X(r, 8) takes a flat reference plane 
containing a circular pore into a curved nearly flat membrane 
with a hole of nearly circular projected area. 



Substituting Eas. lfT^ and ifTH) into Eq.Q) and keeping 
terms up to quadratic order in ff and h, we find 



n = a^{Ll-rl) 



(15) 



- / drdOr 



drdO ra 



V -ff+V ■r]2 . 



The last integral in the above Hamiltonian can be con- 
verted into a line integral by the application of the diver- 
gence theorem. Using Eqs.©, CDl, and ifT^ we arrive, 
after some calculation, to the following form: 



n 



(nr{L: 



(16) 



— / drdO r 
2 



aiVhf + K{\'^hy 



d9 [f-r^{ro,e)]^-aro / d9f-7j{r„,0) 



An important feature of Ea. (|16|l is that, to quadratic or- 
der, the dependence of the Hamiltonian on h{r, 9) and 
f]{r,6) decouples completely. Also, notice that Ea. ljltil) 
depends only on the boundary value of ff. This is 
a consequence of the fact that away from the bound- 
ary, rf is merely a transformation of coordinates under 
which the Helfrich Hamiltonian should be invariant. We 
will henceforth use the scalar function r]{9) to denote 
the the boundary values of the mapping f]{r,9), i.e., 
r]{9) = f ■ fi{rQ,9). Using BC ®, we find out that the 
last two terms of Ea. ljKill cancel each other. This leaves 



us with Ti. ~ cnr{L 



Ti-h = 



drdd r 



Tih , where 



a{\/h:f + n{\/'^hy 



(17) 



The Laplacian in the height-dependent Hamiltonian 
(117(1 requires that we have two BCs on each boundary. 
On the outer boundary (r = Lp) we impose the BCs 



h{Lp) = 0, and, V^h{Lp) = 0. 



(18) 



The first BC corresponds to a membrane which is at- 
tached to a static frame on its external perimeter. The 
second is obtained by considering the discrete version of 
the Helfrich surface Hamiltonian and requiring that in 
the continuum limit, the same equation describes the mo- 
tion of the boundary and bulk elements. The BCs on the 
inner boundary (r — tq) are quite similar: 



h{rQ,9)=H{9) and, V^h{ro)^0, 



(19) 



with the only difference that the height is set to H{9) 
rather than vanishes. The vector 



Y{9) = X{ro,9) = [ro + v{0)] f + H{9)z, 



(20) 



depicts the locus of the pore boundary in the 3D em- 
bedding space. Note that in the case of a membrane 
inclusion of radius ro, the second BC in Ea. (|19() should 
be replaced by -n ■ Vh{ro,9) = dh{rQ,9)/dr = H'{9) 
where the contact slope, H' , depends on the geometry 
and the tilt angle of the inclusion _37J . 

We proceed in analyzing the area term by writing the 
height function as h — ho + f, where ho is the extremum 
of Hamiltonian (|17|l . i.e., 

- crV^/io + kV/io = 0, 

subject to the BCs that 

ho{Lp)=0, ^^ho{Lp) = 0, 



(21) 



hoiro,9)=Hi9), V2/io(^o) = 0. 



(22) 



Eas. H19|l and ((22(1 imply that the function /, which de- 
picts the fluctuations around the equilibrium profile ho 
satisfies 



f{Lp)=0, VV(ip)=0, 
/(ro)=0, v2/(^o) = 0. 

Hamiltonian ((17(1 can be thus written as 



(23) 



nh {ho + 1) = 



drd0 r I — 
2 



(j{Vhof + n{\/'^hof +(24) 



[n^ho ■ V/ + «V2/ioVV] + \[c^ (yff + « (VV)' 



For the cross term (second term in TLh) we obtain, upon 
integration by parts, 

/ drd9 r [crV/io • V/ + kV^/iqVV] = 



f drde r [-crV^/io + kV/io] / - 



df 



27r 

dr 

2-iT 



d [kV^/iq - cr/io] , 

dO J. 7, 

or 



(25) 



where the last two integrals in the above equation are 
performed on the boundaries of the system. The bound- 
ary terms in Ea. lf^ vanish due to the BCs (^3 and (|23J|, 
while the bulk term vanishes due to Ea. (|21l) . 

Without the cross term in Ea. (|24() . the height- 
dependent Hamiltonian takes the simple form 7i^(/io + 
/) = 'Hh{ho) + 'Hh{f), where the energies associated with 
ho (the equilibrium term) and / (fluctuation term) com- 
pletely decouple. Integrating both terms by parts twice, 
we find expressions similar to Ea. H25l) . where ho is re- 
placed by / (in the fluctuation term) or vice versa (equi- 
librium term) . In the former, the boundary terms vanish 
due to Ea. (|23|l and we are left with 



-Hhif) 



drdO r- 



1 



a(V/)' + ^(vV)' 



(26) 



In the latter, the bulk term is eliminated by virtue of 
Ea. (|21|l . Considering the BCs on ho (|22|l . one can easily 
find that 



■Hhiho) 



2. g 

deroH{9) — 



dr 



(27) 



In a manner similar to the last two terms in Eq. Hl()|) , the 
last result demonstrates that contributions to the Hamil- 
tonian due to the pore can only appear through boundary 
(line) integrals. 



C. Pore line tension 

An additional contribution to the free energy is due 
to the line tension of the pore, which arises from the 
curvature and packing of the lipid molecules at the pore 
boundary (fig. Ql. The pore's shape is depicted by the 
curve Y(0) (|20|l . and the line tension energy is given by 



(.27r 

r / de 

Jo 

r r^ 

^^oJo 



{dY{0)/d0y ~ r 27rro + T dO 77(6*) 



2ir 



dO 



{dH{e)/def + {dr]{9)/dey 



(28) 



where T is the line tension coefficient, i.e., the edge energy 
per unit length. 



D. The full Hamiltonian 

Collecting expressions ((THJ, ^^, (EH), and ((^HJ, we 
find that the Helfrich Hamiltonian (excluding the Gaus- 
sian curvature term) can be written as the sum of two 



terms: 



H = Hs +^1- 



(29) 



The first term is the surface Hamiltonian associated with 
the membrane fluctuations: 



Hs = CFT^iLp 



drdO r 



<y{^lf 



i^'fY 



(30) 

The second is the inner boundary term, consisting of the 
various contributions to the energy due to the pore 



Hi = T2TTro 

1 r2iT 



de 







--^(^)' 



(31) 



+ 



r (dH{9)Y 

ro\ de J 



roH{e) 



d {nS/^ho — crho) 



dr 



III. SURFACE TENSION 

In calculating the contribution of the surface tension 
term in Hamiltonian (|30|l to the free energy, we follow the 
procedure described in our previous manuscript on mem- 
brane inclusions 37J . We start by expanding the function 
/ in a series of eigenfunctions fmni'r) of the operator 
C = -aV^ + nV^: /(r,0) = E™,n^™,n/™-(r)e™f 
The functions fm,n{i') can be written as the linear com- 
bination of the Bessel functions, Jm{r) and Ym{r), of the 
first and second kinds of order to, and the modified Bessel 
functions of the first and second kinds of order to, Kmij") 
and Im{r): 

jmAr) = AJ„(Ar'"r)-Hi?i;„(Ar'"r) 
+ CKrn{XT''r)+DIrniXT''r), 

where the Xi {i — 1,2) are the positive solutions 
of (-1)*+V(A™'")2 -I- k(A™'")'* = ij,m,n, and ii^.n is 
the eigenvalue corresponding to the function fm,n{f)'- 

Applying the BCs it^ at tq and Lp, we derive the 
eigenvalue equation 

Jm(A'"'"ro)Fm(A7^"ip) - J„i(A7'"-Lp)i^m(A"'"ro) = 0. 

(32) 
Although this eigenvalue equation is different from the 
one we had in Ref. 37], the asymptotic behavior of the 
eigenvalues is the same. In the long wavelength limit, 
A™'"ro ^ |to|, Eq.(|22Jl reduces to the eigenvalue equation 
in the absence of pores: 



Jm\Xi 



'Lp) = 0. 



(33) 



This is a manifestation of the fact that modes with char- 
acteristic lengths much larger than the pore radius are 
hardly perturbed by its presence. In the opposite limit, 
A™'"?'o ^ |'7i|, we find that the difference between two 



consecutive eigenvalues saturates to Aj^ — A^^ — 

iT/{Lp — vq), which is a factor Lp/{Lp — vq) larger than 
in the case with no pore. The physical interpretation 
of this result is that the pore acts like a hard wall for 
modes with characteristic lengths much smaller than its 
radius, reducing the effective linear size of the membrane 
for these modes to Lp — rg. 

The insertion free energy associated with the surface 
Hamiltonian, defined as AFs{ro) = Fs{ro) — Fs{0), can 
be expressed by the following sum [33 



AF,{ro) 



-naorl 



(34) 



El- 



/ \m,n \9 I / \m,n \a 



LI 



where the sum runs over the modes n = 0, 1, ... , y/No, 
and, 171 = — \/No, . . . , ^/Nq so that the total number of 
modes 2A^o is proportional to the number of molecules 
forming the membrane N, while A™;q> are the corre- 
sponding solutions of the eigenvalue equation in the ab- 
sence of pores (p^ . As in j^^j Ea. |Mjl was derived 
by expanding the surface tension to quadratic order in 
ro/Lp, which is valid for small pores tq <C Lp. Ana- 
lytical approximation of this expression is obtained by 
assuming (based on our discussion of the asymptotic be- 
havior of the eigenvalues Aj,"'") that eigenvalues such that 
A™'"ro < a\m\ (long wavelength) are not affected by the 
pore, whereas modes with A™'"ro > a\m\ (short wave- 
length) grow by a factor Lp/{Lp~ro). The dimensionless 
constant a is of the order of unity and its value will be 
determined later by exact numerical calculation of AFg . 
Using this "step-function" approximation for the eigen- 
values A™'", and evaluating the sum in Ea. (|34|l as an inte- 
gral, we obtain the simple result (correct up to quadratic 
order in ro) that 



AF, ^ Fs{ro) ~ Fs{0) = 




9 ksTrl 





In 



= -vrro(cr + Act) = -Uriahs, 



(35) 



where S, = \J nja, and Iq = Lp/^/N^ is a microscopic 
length cutoff which is of the order of the bilayer thick- 
ness. From the above equation, we identify the thermal 
correction to the surface tension as 



AtT = 



keT 
nalg 



-^1 lu 



if' +' 



(36) 



In order to test the accuracy of expression (|35|l . we 
have numerically solved the eigenvalue equation (|32|l and 
used the solutions to evaluate the sum in Ea. H34|l . Nu- 
merical values of AFs{ro) (for k = IO/cbT and various 
values of a) are shown in fig. |3] (a)-(b). They have been 
extracted by extrapolating the numerical results obtained 
for several values of 750 < A"o < 2000 to the thermody- 
namic limit Nq — > oo. In the inset to fig. |3| (a), the 
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FIG. 3: The surface free energy AFs as a function of the 
inclusion's radius for n — lOfcsT and various values of a. 
The inset to graph (a): a log- log plot of the numerical results 
for a — 0. The slope of the straight dotted line is 2. 



results for a — are replotted on a logarithmic scale, 
showing that our prediction of a quadratic relation be- 
tween AFs and tq is attained only for large (macroscopic) 
pores with rg > IOOZq (the slope of the straight dotted 
line is 2). The discrepancy between the numerical val- 
ues of AFs and Ea. (|35|l in the small rg regime is due 
to the significant contribution to the free energy of the 
crossover modes A™'"ro ^ 1 which is poorly calculated 
by the "step function" approximation. The solid curves 
in fig.|Sl(a)-(b) depict our analytical expression (|^ for 
AFs , with a determined by fitting the results for large tq 
to Eq. (|35|) . The value of a shows a slight dependence on 
the surface tension varying from 1.60 for a — to 1.75 
for ^ — a/k/ct — ^Iq/tt. 

Our numerical and analytical results suggest that 
Act < 0, making the effective surface tension smaller than 
the bare surface tension. Of particular interest is the 
fact, demonstrated in fig.|Sl(a), that for weakly stretched 
membranes (large ^) the effective tension may be nega- 
tive. In such a case the effective surface tension would act 



to prevent, rather than facihtate, the opening of a pore. 
For strongly stretched membranes (small £,) the domi- 
nant contribution to a^s is of the bare surface tension. 
In this regime, the surface tension part of the free energy 
is well approximated by the second term of Eq.(^, i.e., 
SFs ~ -anr^. 



IV. LINE TENSION 

In order to calculate the contribution of the line tension 
term in Hamiltonian H31|l to the free energy, one needs to 
trace out the fields rj and /ig which are decoupled from 
each other. Introducing the Fourier transform of the field 
V, 



vio) 



MB 



m=-Ni 



(37) 



where ddB is the de-Broglie thermal wavelength. Making 
the particular choice of 770 that satisfies BC ©, 
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MB 
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\Vr, 






(38) 



the corresponding Hamiltonian takes the form 



that (for a > 0) ho can be written by the following mode 
representation 



ho{r,9) 



MB 



-I- ^n— _ AT, 



irnO I r 



m=-N 



In (r/rp) 
In (ro/Lp) ' 

(43) 
Notice that V^ft-o — everywhere and not only at the 
boundaries. Substituting expression H43|l in Ea. H31|) . one 
arrives to the following Hamiltonian 
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tq 



(44) 



Assuming that a <C F/rp (the weak stretching regime), 
it is easy to conclude that the resulting contribution to 
the free energy Fl^ ~ F^ , and thus 



(45) 



= 2^ro(r + AF) = 27rroreff. 

We thus identify the thermal correction to the line ten- 
sion of the pore 



Fi - 


= 2^ror + F/' + F;'' 
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(39) 



The |?7i| = 1 modes are trivial translation modes which 
do not contribute to the energy of the pore. The number 
of m > 2 modes is equal to the number of microscopic 
degrees of freedom, namely the number of molecules on 
the rim of the pore. Since this number is proportional to 
the perimeter of the pore, we can write 



A^i 



(40) 



where 6 is a numerical factor of the order of unity. Trac- 
ing over the variables fim is straightforward, giving the 
free energy 



FF 



knT 



El- 



r/i > 1 



dlB 



- 1) r/ro 



keTNi 



(41) 



If the number of modes is large A'^i ^ 1 (i.e., lo <C tq) 
then the sum in the above expression can evaluated as 
an integral, giving 



F.^ ~ 27rro 



bkBT 



2 [ksTlo 



-1 



(42) 



The contribution of the field Hq to the line tension 
free energy is also tractable. From the partial differen- 
tial equation H21|) and the BCs H22|) . it is easy to show 



AF 



bkeT 
ttIq 






(46) 



For phospholipid bilayers at room temperature AF is neg- 
ative and is typically in the range of 10~^ — 10~^ dyn. 



V. DISCUSSION AND SUMMARY 

The opening of a membrane pore has been tradition- 
ally regarded as an energetically-driven process. Accord- 
ing to this view, the surface and line tensions are the 
forces driving, respectively, the opening and closure of 
pores. The balance between these opposing forces cre- 
ates a nucleation barrier for the formation of long-lived 
pores, and requires the opening of a sufficiently large hole 
at the initial stage. 

In previous studies, the role of thermal fiuctuations 
has been limited to facilitating the opening of a nucle- 
ation pore. The critical pore size ^ and the height of 
the barrier @ have been determined from Eq. ^ for the 
pore energy. However, at nonzero temperature an en- 
tropic part must be added to the Eq.Q. To fill the gap 
in the literature on the subject, we have calculated the 
thermal contributions to the pore free energy associated 
with (a) the shape of the boundary of the hole, and (b) 
the fluctuation spectrum of the membrane. Our study 
suggests that the pore free energy may be expressed by 
an equation similar to (^ 



F = Feff27rro - dcSTrrg, 



(47) 



8 



in which the bare surface and hne tensions are replaced 
by effective (renormaUzed) values. Typically, we find that 
Fcff < r and cTcff < cr, reflecting two opposite tenden- 
cies. The decrease in the line tension reduces membrane 
stability against pore formation. It reflects the larger 
configuration space available to a membrane with a hole 
present. The decrease in the surface tension, on the other 
hand, makes the formation of pores harder in comparison 
to the zero-temperature case. This effect originates from 
the change in the spectrum of membrane fluctuation oc- 
curring upon the opening of the pore and the resulting 
increase in bending energy. 

We can identify a number of different regimes of pore 
stability. For tense membranes with positive effective line 
tension, we find the standard regime of classical nucle- 
ation theory. From Eq. H47|) , we can identify the stability 
criteria for the growth of large pores to be tq > Fcft /coff . 
For membranes with low surface tension, the effective 
surface tension will be negative. In this regime, pores 
will increase the free energy for all radii and one should 
not expect the formation of pores spontaneously as long 
as the effective line tension is positive. This regime is 
quite unlike classical nucleation theory, where a nucle- 
ation barrier for pore formation always exists. 

In the theory of thermally activated poration, the nu- 
cleation rate of critical pores depends strongly on the the 
free energy barrier AF = 7tTcs/<JcS, as exp{~ AF / ksT) . 
The height of the nucleation barrier decreases as one ap- 
proaches the temperature at which the effective line ten- 
sion, Foff H45|l . vanishes. Above this temperature the bar- 
rier disappears and the formation of pores occurs spon- 
taneously. The growth of the pore will not be stopped 
as long as the surface tension remains positive. If the 
rate at which the surface tension is relaxed is too small, 
the membrane will rupture. However, we have demon- 
strated in section ITTll that membrane fluctuations renor- 
malize the surface tension, making the effective tension 
negative when the applied tension is small. Therefore, we 
may have a situation where both Foff and (t^s are nega- 
tive. In such a case, the free energy attains a minimum, 
not a maximum, for ro = Foff/tToff, and a stable pore of 
that radius will be spontaneously formed. The radius of 
such a pore can be easily varied by several orders of mag- 
nitude by varying the applied tension and thereby tuning 
the value of a^s- 



The pore size can be also varied by manipulating the 
magnitude of the line tension, e.g., by the addition of 
colipids that modify the ability of the lipids to pack at 
the edge region. Experiments in bilayer lipids to which 
lysoPC and cholesterol were added demonstrated that 
these molecules affect the line tension in opposite ways, 
with the former decreasing the line tension |39l] and the 
latter increasing it 23, 25]. Typical values of the line 
tension found experimentally are in the range of our es- 
timate of the thermal correction AF or somewhat larger. 
This demonstrates the significance of thermal effects in 
determining the stability of membranes against pore for- 
mation. Reexamining previous studies which do not in- 
clude temperature-dependent corrections is of particular 
importance because the experimental determination of 
the line tension T^s is not direct. It is based on a suit- 
able model relating the line tension to other measurable 
quantities such as the surface tension, pore radius, etc. 

To better understand the effect of thermal fluctuations 
on pore formation, it is necessary to go beyond the sec- 
ond order expansion of the Helfrich Hamiltonian in h and 
■q. In the Foff < regime, the system gains free energy 
by having pores with large perimeter, which means that 
holes with shapes that strongly deviate from circular will 
be highly favorable [33 • By restricting our analysis to 
small values of the field rj (and hence to quasi-circular 
pores only), we probably underestimate the magnitude 
of the thermal correction to the line tension. Moreover, 
for negative values of the effective line tension the forma- 
tion of many pores is also likely to occur. The collective 
behavior of these pores and the (membrane-mediated) 
interactions between them w| arc not well understood. 
More insight on the pores' architecture and their dynami- 
cal evolution can be gained by molecular- level studies and 
computer simulations. 
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